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晶体 相 场 法 研究 Kirkendall 效 应 诱发 的 
相 珊 空洞 的 形成 和 演变 - 
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摘 要 采用 二 元 晶体 相 场 模型 研究 了 Kirkendall 效应 诱发 的 相 界 空洞 的 形成 及 扩展 过 程 . 模拟 结果 表明 : 对 于 取 
较 小 的 相 界 , 空洞 向 原子 迁移 率 高 的 一 侧 (a 相 ) 方向 移动 , 空洞 的 形状 由 最 初 的 平行 四 边 形 向 六 边 形 演化 , 空洞 周围 原 
子 淹没 速率 大 于 产生 速率 , 从 而 造成 空洞 扩大 , 空洞 扩展 过 程 中 伴 有 相 界 的 移动 及 相 的 长 大 和 缩小 . 对 于 取向 差 角 较 
的 相 界 , 空洞 还 会 沿 相 界 方向 扩展 , 使 得 空洞 连通 , 将 相 界 割裂 开 , 割裂 后 的 两 侧 相 界 呈现 锯齿 状 . 扩散 过 程 中 , 体系 的 
自由 能 逐渐 降低 . 对 于 取向 差 角 较 小 的 相 界 , 原子 迁移 率 差 值 增 大 , 自由 能 下 降 无 明显 差异 . 对 于 取向 差 角 较 大 的 相 界 ， 
原子 迁移 率 相 差 越 大 , 自由 能 下 降 速 率 越 快 . 随 着 相 界 取向 差 角 的 增 大 , 自由 能 的 下 降 速率 逐渐 增 大 . 相 界 空洞 的 模拟 
结果 与 实验 观察 一 致 . 
关键 词 二 元 晶体 相 场 模型 , Kirkendall 效 应 , 相 界 , 空位 , Kirkendall 空洞 
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ABSTRACT The mechanical properties of materials are related to the integrity of interfaces (phase and grain 
boundaries). For substitutional alloys, the Kirkendall voids tend to form more easily at the phase boundary or grain 
boundary when the atomic mobilities of different species are unequal, which will degrade the bounding quality of 
interfaces. So far, there have many experimental studies on the evolution of Kirkendall voids and the formation 
mechanism. However, allowing for the fast process of the Kirkendall voids from formation to evolution, it is hard 
to capture such process in real experimental conditionals. So the formation and evolution mechanism of the Kirken- 
dall void need to be studied. A binary phase field crystal model was used to simulate the process of void formation 
and expansion at phase boundaries induced by the Kirkendall effect. Simulated results show that for the low misori- 
entation phase boundary (PB), the void moves toward the side with large atomic mobility (& phase) and the void 


shape evolves from the initial parallelogram to hexagon. The atomic annihilation rate around a void is faster than 
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that of growth rate, which results in void expansion. The PB migration, phase growth and shrinkage can also be ob- 


served in void expansion. For the large misorientation PB, voids can also expand along the PB direction, resulting 


in the connection of voids, therefore, the PB is separated and presents zigzag shape. In the interdiffusion system, 


the free energy decreases. The descending speed of the free energy 1s almost equal for the low misorientation PB 


while is increasing for the large misorientation PB when the atomic mobility difference becomes larger. The de- 


scending Speed of the free energy is proportional to PB misorientations. The PB void predicted from our computer 


simulation is consistent with the experiment observation. 


KEY WORDS binary phase field crystal, Kirkendall effect, phase boundary, vacancy, Kirkendall void 


材料 的 力学 性 能 与 界面 ( 晶 界 、 相 界 等 ) 完 整 性 
有 关 , 影响 界面 质量 的 因素 有 很 多 , 如 杂质 元 素 偏 
聚 、 晶 界 预 熔 等 . 除 此 之 外 , 对 于 置换 型 合金 , 晶 界 
或 相 界 处 往往 会 出 现 空洞 , 这 些 空 洞 称 为 Kirkendall 
空洞 中 会 显著 影响 界面 的 质量 , 降低 材料 的 力学 性 
能 . 在 电子 封装 领域 , Kirkendall 空洞 的 危害 尤为 突 
出 . 随 着 电子 元 器 件 向 小 型 化 及 多 功能 化 发 展 , 封 
装 密度 不 断 提 高 , 要 求 焊料 连接 处 的 尺寸 不 断 减 
小 . 电子 的 迁移 及 元 器 件 长 时 间 的 高 温 工 作 状 态 使 
得 焊料 在 连接 处 产生 Kirkendall 空洞 ”大 大 缩短 
了 电子 元 器 件 的 寿命 . 这 种 Kirkendall 空洞 是 焊料 原 
子 与 被 焊 金 属 原子 的 非 等 量 扩散 造成 的 . Kim 等 % 指 
出 , 在 Cu/CusSn 扩散 偶 的 互 扩 散 过 程 中 , Cu 原子 是 
CuwsSn 相 及 CucSns 相 中 的 主要 扩散 原子 . 相 界 处 , 由 
于 Cu 原子 的 迁移 能 力 大 于 Sn 原子 的 迁移 能 力 , 发 
生 不 等 量 原子 扩散 , Sn 原子 无 法 及 时 弥补 由 于 Cu 
原子 扩散 产生 的 空缺 位 置 . 这 样 , 一 个 原子 尺度 的 
空位 便 被 保留 下 来 , 随 着 互 扩散 过 程 的 进行 , 空位 
发 生 集聚 (多 发 生 于 位 错 等 缺陷 处 ), 形成 空洞 . 迄今 
为 止 , Kirkendall 空洞 的 形成 及 演化 机 制 已 经 有 了 许 
多 实验 方面 的 研究 59. Miiller 等 "提出 了 Kirkendall 
空洞 形成 及 扩展 的 连续 模型 . Nakahara 和 McCoy™ 
研究 了 注 层 扩散 偶 Kirkendall 空洞 的 形成 过 程 , 发 
现 空洞 与 扩散 偶 的 厚度 .形状 及 唱 粒 尺寸 有 关 . 
Zeng 等 在 研究 SnPb 焊料 应 用 于 电子 封装 领域 时 ， 
发 现在 高 温 与 承受 应 力 的 情况 下 , 在 CuSn 与 Cu 界 
面 附近 观察 到 了 大 量 的 Kirkendall 空洞 , 降低 了 焊 
接 的 牢固 性 . Son 等 吧 对 Sn:Ag/Cu 扩 散 偶 进 行 热 老 
化 实验 (恒温 加 热 的 方式 来 单独 研究 热 场 对 接头 可 
靠 性 的 影响 ) 时 , 发 现 随 着 热 老 化 时 间 的 延长 , 相 界 
处 的 Kirkendall 空洞 会 发 生 连 接 . Xu 等 趾 对 出 现 
Kirkendall 空洞 的 扩散 偶 进 行 冲击 实验 , 结果 表明 
Kirkendall 空洞 会 降低 扩散 侦 的 冲击 强度 . 尽管 上 述 
研究 在 一 定 程 度 上 对 Kirkendall 空洞 的 形成 机 制 做 
出 了 解释 , 但 是 由 于 Kirkendall 空洞 的 形成 到 长 大 
是 一 个 极为 快速 的 过 程 , 因此 很 难 揭 示 出 Kirkendall 


空洞 的 形成 及 演化 规律 . 除 此 之 外 , Kirkendall 空洞 
的 出 现 对 于 实验 条 件 的 要 求 也 较为 苛刻 , 实验 过 程 
中 很 难 排除 外 在 因素 (如 杂质 原子 、 退 火 温度 及 扩散 
偶 的 类 型 等 ) 对 Kirkendall 空洞 形成 及 演化 的 影响 . 
实验 的 周期 较 长 , 花费 较 大 , 不 可 控 因 素 多 等 都 使 
得 对 于 Kirkendall 空洞 的 实验 研究 较 难 开展 . 

近 几 年 随 着 计算 机 技术 的 发 展 , 数值 模拟 技术 
得 到 了 长 足 的 发 展 , 可 以 模拟 材料 微 结构 之 间 界 面 
( 晶 界 、 相 界 等 ) 的 形 貌 , 而 这 将 有 助 于 揭示 相 界 处 
Kirkendall 空洞 的 形成 及 演化 规律 . 然而 , 迄今 为 目 
对 于 Kirkendall 空洞 方面 的 计算 机 模拟 研究 还 很 缺 
乏 . Wu 等 中 利用 相 场 方法 模拟 了 二 元 合金 中 纳米 
空洞 的 形成 过 程 . Delogu05 采 用 Monte Carlo 方法 研 
究 了 纳米 空洞 的 稳定 性 问题 . 杨扬 中" 采用 第 一 性 原 
里 的 方法 研究 了 CuSn 相 及 CusSn: 相 中 的 空位 扩散 
性 质 . 

尽管 上 述 方法 可 以 在 一 定 层面 上 反映 Kirkend- 
all 空洞 的 相关 问题 ,但 是 上 述 研究 体系 均 无 法 综合 
考虑 弹 塑 性 、 各 向 异性 、 晶 粒 取 问 等 诸多 物理 特性 
对 扩散 体系 的 影响 , 有 关 Kirkendall 效应 诱发 的 相 
界 空洞 的 模拟 研究 还 很 缺乏 . 

晶体 相 场 (phase field crystal, PFC) 模 型 点 "9 是 在 
密度 泛 函 理论 (DFT) 基 础 上 发 展 起 来 的 , 耦合 了 包括 
弹 塑 忻 、 各 向 异性 ` 唱 粒 取 向 等 诸多 物理 特性 , 可 以 
在 更 大 的 时 间 尺 度 (扩散 时 间 尺 度 ) 和 更 小 的 空间 尺 
度 (原子 空间 尺度 ) 上 模拟 材料 的 相关 行为 , 已 经 被 成 
功用 来 模拟 晶 界 运动 ”以 及 材料 裂纹 扩展 中 的 相 
关 缺 陷 行 为 . 本 工作 采用 三 元 晶体 相 场 模型 ”1 研 
究 了 扩散 过 程 中 空洞 的 形成 与 扩展 过 程 , 重点 分 析 
了 其 与 Kirkendall 效应 的 关系 . 

1 模拟 方法 
1.1 二 元 晶体 相 场 模 型 

基于 密度 泛 函 理论 的 二 元 晶体 相 场 模型 由 El- 
der 等 中 提出 , 该 模型 用 2 个 场 变 量 来 表达 自由 能 函 
数 , 即 原子 密度 场 n= (p+ps -Pp )/p' , 原子 浓度 场 
w=(ps -ps)/P (其 中 pw， ps 和 p 分 别 为 A 原子 、B 
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原子 与 参考 液 相 的 原子 密度 ). 在 上 述 条 件 下 , 忽略 
原子 尺度 差异 , 置换 型 二 元 唱 体 相 场 模型 的 自由 外 
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式 中 , MA 和 Ms 分 别 为 A, B 2 种 原子 的 原子 迁移 
率 . 基于 时 间 采 用 的 是 无 量 纲 形式 , 为 了 简化 计算 ， 


无 量 纲 形式 所 表述 如 下 : 
加 nA _t,3, 0 4 
F= Jar|34m 3 二 


+ 号 + 全 (GD 
式 中 ，F 表征 无 量 纲 的 自由 能， A"= 
AB,+ BWP+B;(1+) ， AB。=B,-B; 表示 无 量 纲 
的 过 冷 度 ，B 表征 无 量 纲 的 液 相 体积 模 量 ，B; 正比 
于 固 相 晶体 的 弹性 常数 ，B, 是 无 量 纲 的 液 相 体 积 模 
量 Taylor 展 开 式 中 的 第 二 项 系数 , 决定 了 液 相 体积 
模 量 随 原 子 浓度 场 的 变化 情况 ; K 表征 浓度 场 的 梯 
度 能 系数 ; w 决定 了 液 固 两 相 之 间 相 互 转换 的 趋 
势 ; u ,与 w@ 共 同 决定 了 相 分 离 的 尺度 ; 1, vv 与 y 
是 模型 参数 ; 在 保证 不 失掉 二 元 晶体 相 场 模型 基本 
物理 特性 的 前 提 下 , 为 了 简化 计算 , 在 本 工作 中 参 


数 y=0. 
A, B 2 种 原子 的 密度 场 演 化 方程 中 如 下 : 
On _ 2 OF _ 
本 


MV 4 -tn ton+ [+ By + uy Kv2y] (2) 
One _ 2 OF _ 
Msvi| An -tn +vn — [eo + Bly -w+ Kvy| (3) 
根据 原子 密度 与 浓度 场 的 关系 n= 
(n+y)/2=pa/p 及 ns=(n-y)/2= pa/pi， 可 以 得 
到 n=n、+ns 及 小 =n, -ns. 经 过 推导 可 以 得 到 原子 
密度 场 n 及 浓度 场 y 的 动力 学 演化 方程 : 
on _ an tns) _ ona 下 Ons _ 
Of 0 ot 0 
2 OF 26F _ 
MV Bn Br 


(M, 十 Ms)V (An —in + vn’) 十 
(WM, -Ms) Vl(w + Bs)y + up|+ 
(Ms -MA)V KVY (4) 


(Ms Ms)V (A'n —in + wn’) 十 
(0 +HMnvi(w +B)y+ wa _ 
(Ms + MA)V KV (5) 


设 定 A 原子 迁移 率 为 1 ( 即 Mo=1, 且 Mr>Ms). 对 于 
上 述 动力 学 方程 时 间 方 向 的 求解 采用 Euler 方法 ， 
空间 方向 的 求解 采用 半 隐 式 Fourier 谱 方 法 后 
离散 后 使 体系 进行 演化 , 得 到 原子 密度 及 原子 浓 
度 信息 . 
1.2 模拟 参数 及 初始 设置 

本 工作 中 采用 的 二 元 晶体 相 场 模型 参数 的 
选取 参考 Elder 等 所 的 工作 . 模型 参数 选 定 为 : 
有 =0.7, B, =-1.8,， Bi; =1, 大 0.6, v=1, K=4, @=1, u 
4. 初始 的 原子 密度 设置 为 : 固 相 n(x,y,t1=0)= 
元 十 Heos(24,7)/2 = cos(g,x)eos(q,9)) = -0.1501 + 
Becos(24,7)/2 一 cos(gszjcos( 9， 7 , 液 相 n(x,y,t=0)= 
而 =-0.2519, 其 中 4.=27 9， = a 为 原子 蝇 
格 常数 , ,与 元 分 别 为 固 相 与 液 相 中 初始 原子 密度 
平均 值 . 在 过 冷 液 相 中 设置 一 个 晶 粒 , 在 演化 公式 
(4) 和 (5) 的 驱动 下 , 计算 体系 进行 弛 豫 , 体系 达到 液 - 
固 平 衡 状态 , 从 而 得 到 平衡 后 液 相 与 固 相 的 平均 原 
子 密度 分 别 为 : 元 =-0.2517, 元 =-0.1503 . 平衡 唱 
格 常数 为 : a =6.927. 计算 区 域 为 [x=1024Ax 及 
Ly=1024Ay , 采用 周期 性 的 边界 条 件 . 空间 步 长 及 
时 间 步 长 分 别 为 Ax=Ay=a./7 及 At=0.1. 初始 浓 
度 场 的 设置 如 图 1 所 示 . 


~ 


NA 


NA 


1 Kirkendall 效 应 诱发 相 界 空洞 的 初始 设置 示意 图 
Fig.1 Initial setup schematic of phase boundary void in- 
duced by the Kirkendall effect (x and y indicate the 


two directions of the Simulated area. w, PB, J,, Je 
and 和 indicate the atomic concentration field, the 
position of the phase boundary, the A atomic flux 
from & phase to bP phase, the B atomic flux from 
phase to w phase and the net vacancy flux from £6 


phase to phase, respectively) 
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2 结果 及 讨论 

在 相 界 空洞 形成 及 扩展 的 模拟 中 , 计算 区 域 由 
2 个 不 同 取向 差 角 的 相 区 组 成 , 并 且 上 部 a 相 设 置 
初始 浓度 y=0.15 ( 富 A 原子 区 域 ), 下 部 8 相 设置 初 
始 浓度 =-0.15 ( 富 B 原子 区 域 ). 

模拟 开始 阶段 , 首先 设置 2 类 原子 迁移 率 为 
Ms=Me=1, 并 在 不 同 取 向 的 两 相 之 间 设 置 很 薄 的 过 
冷 液 相 , 使 之 发 生 凝 固 从 而 在 计算 区 域 中 间 位 置 得 
到 相 界 , 如 图 2 和 3 中 黑色 虚线 (o ) 所 示 . 随后 , 调整 
原子 迁移 率 为 M,=100M, 使 体系 发 生 互 扩散 . 为 了 
排除 周期 性 边界 条 件 产生 的 另 一 相 界 的 影响 , 仅 选 
取 计 算 区 域 的 一 部 分 ( Lx=100Ax, Ly=512Ay ) 进 行 
分 析 . 
2.1 取向 差 角 较 小 的 相 界 空洞 

图 2 为 取向 差 角 较 小 的 相 界 (2°) 空洞 的 生成 
及 扩展 过 程 . 从 图 2a 可 以 看 出 , 当 演 化 时 间 
T=1x10At 时 , 靠近 w 相 一 侧 (高 原子 迁移 率 ) 有 一 
空洞 形成 , 形状 呈 平 行 四 边 形状 , 在 随后 的 扩展 过 
程 中 (图 2b~e), 空洞 向 a 相 一 侧 发 生 移动 , 并 且 不 断 
扩大 , 形状 由 平行 四 边 形 逐 渐 演变 为 六 边 形 . 空洞 
移动 过 程 中 , 相 界 也 向 a 相 一 侧 发 生 移 动 , 即 6 相 长 
大 , w 相 缩小 . 

Kirkendall 效应 发 生 在 置换 型 合金 的 原子 扩散 
过 程 中 , 原子 通过 与 空位 发 生 位置 交 换 , 实现 扩散 
过 程 . 由 于 不 同 原 子 的 原子 迁移 率 存在 差异 , 空位 
首先 倾向 于 与 迁移 率 大 的 原子 发 生 位 置 交 换 , 因此 
在 整个 扩散 过 程 中 存在 净 原 子 流 和 净空 位 流 . 基于 
空位 机 制 , 体系 扩散 过 程 中 , 始终 满足 = 及 .= 
wu (其 中 为 A 原子 由 a 相向 6 相 的 扩散 通 量 ; .NA 


(a) (b) (9) (d) (e) 
00 
9 一 一 一 
y | 
b 


2 取向 差 角 较 小 的 相 界 (0=2°) 空 洞 的 形成 及 演化 


Fig.2 Void formation and evolution diagrams at the low 


一 


misorientation (=2? ) phase boundary at evolution 
time 7=1 x 10°At (a), T=2x 10°At (b), 7=3 x 10°At (¢), 
T=4 x 10°At (d) and 7=6 x 10°At (e) (o indicates the 
initial position of the phase boundary) 


为 对 应 于 A 原子 的 扩散 通 量 并 与 其 方向 相反 的 空 
位 扩散 通 量 ; .为 B 原 子 由 6 相向 a 相 的 扩散 通 量 ; 
Ju 为 对 应 于 B 原子 的 扩散 通 量 并 与 其 方向 相反 毕 
空位 扩散 通 量 ). 

本 部 分 中 , 取向 差 角 较 小 的 相 界 可 以 视 为 由 间 
距 较 大 的 一 组 位 错 核 组 成 , 从 图 2a 可 以 看 出 , 在 选 
定 的 分 析 区 域 仅 包含 一 个 位 错 核 , 该 位 错 核 可 以 作 
为 空位 的 淹没 源 . 在 扩散 过 程 中 , 由 于 空位 首先 倾 
向 于 与 原子 迁移 能 力 强 的 A 原子 发 生 位 置 交 换 , 则 
由 Qa 相向 6 相 一 侧 存 在 净 的 原子 扩散 通 量 A 太 /~ 
太 , 由 相向 a 相 一 侧 存 在 净 的 空位 扩散 通 量 A/= 
Jw-Jve. 空位 扩散 至 空洞 的 上 部 , 与 空洞 上 部 周围 的 
原子 发 生 复 合 反 应 , 原子 位 置 消失 , 该 过 程 使 空洞 
发 生 扩 展 . 原子 扩散 至 空洞 的 下 部 , 并 在 空洞 下 部 
发 生 集聚 , 最 终 使 得 空洞 下 部 周围 析出 原子 , 该 过 
程 使 得 空洞 发 生 淹没 . 由 于 空位 的 潭 没 (复合 反应 ) 
速率 比 原子 集聚 析出 的 速率 要 快 , 故 空洞 在 位 错 核 
的 基础 上 向 a 相 一 侧 移 动 并 不 断 扩大 . 也 正 是 基于 
此 , 6 相 不 断 变 大 而 a 相 逐 渐 缩 小 . 除 此 之 外 , 空洞 
形状 由 平行 四 边 形变 为 六 边 形 是 与 模型 有 关系 的 ， 
该 模型 模拟 的 是 六 角 相 , 故 空洞 尺度 变 大 后 也 倾向 
于 呈 六 边 形 . 
2.2 取向 差 角 较 大 的 相 界 空洞 

3 为 取向 差 角 较 大 的 相 界 (GQ=12°) 空 洞 的 生成 


(a) Nl (b) Ba (c) BR (d) BR (¢) 
[OU 
0O_ 
有 
() BM (0) n) 汪 二 () 0) | 
O_ 


( 
( 


图 3 取向 差 角 较 大 的 相 界 (9-129) 空 洞 的 形成 及 演化 


Fig.3 Void formation and evolution diagrams at the large 
misorientation (=12°) phase boundary at 7=8x10°At 
(a), 7T=6x10Ar (b), T=8x10Af (ce), T=9x 10At (d), 
T=1x10°At (e), T=1.1x10At (), 7=1.3x10°At (g), T= 
1.4x105At (h), 7=1.6x10At Gi) and 7=2x10:At ()) 


202303.00596v1 


chinaXiv 


870 


及 扩展 过 程 由 


图 3a 可 以 看 出 , 当 演 化 时 间 


T=8x10°Az 环 


排 的 2 个 空洞 形成, 形状 星 平行 四 边 形状 , 在 随后 的 
扩展 过 程 中 (图 3a-j), 空洞 向 a 相 一 侧 发 生 移动, 3 
且 不 断 扩大 , 形状 也 由 平行 


, 靠近 w 相 


侧 (高 原子 迁移 率 ) 有 3 


a 
«Iz 


四 边 形 向 六 边 形变 化 . 


除 此 之 外 , 与 取向 差 角 较 小 的 相 界 不 同 , 空洞 之 间 


也 发 生 了 连通 从 图 


T=1.4x10Azs 时 ， 


被 割裂 开 , 两 侧 相 界 呈现 锯齿 状 , 相 界 


3h 可 以 观察 到 在 


由 于 空洞 之 间 的 连通 , 使 得 相 界 


也 会 发 生 类 


似 取向 差 角 较 小 的 相 界 的 移动 , 空洞 扩展 连通 与 相 
界 移动 的 结果 使 得 6 相 长 大 , w 相 缩小 . 

取向 差 角 较 大 的 相 界 空洞 未 连通 之 前 的 空洞 
形状 、 空 洞 运动 及 扩展 类 似 于 取向 差 角 较 小 的 相 界 
的 情况 . 空洞 的 横向 扩展 连通 可 以 解释 为 : 本 部 分 
取向 差 角 较 大 的 相 界 可 以 视 为 由 间距 相对 较 小 ( 达 
到 了 空洞 之 间 相 互 扩散 的 临界 值 ) 的 位 错 核 组 成 , 在 


空洞 的 演化 过 程 中 ,1 


度 差 , 使 得 空洞 之 


于 不 同 空 洞 之 间 存 在 空位 浓 


间 的 原子 与 空位 发 生 复合 反应 (多 


余 空位 


没 ), 从 而 使 空洞 之 间 逐 渐 达 到 连通 . 


2.3 自由 能 与 原子 迁移 率 及 相 开 取向 差 角 的 关系 


为 了 描述 空 ] 


子 迁 移 率 对 体系 自由 能 


的 相 界 取向 差 角 (G=2°, 0-7.6", GQ=12° 及 G24°) 和 原 


影响 , 本 研究 中 选 


同形 成 过 程 中 相 界 取向 差 角 及 原 


了 不 同 


子 迁 移 率 比值 (WM=Ms, Ms,=10Ms, M\=50Ms 及 MA= 


100M) 进行 研究 . 


图 4 反映 的 是 不 同 原子 迁移 率 下 , 取向 差 角 较 
小 的 相 界 (6G=2°) 与 取向 差 角 较 大 的 相 界 (G12°) 扩 散 


过 程 体系 的 自 1 
看 出 ， 


自 


扩散 体系 , 扩散 前 期 , Ms=Ms 对 应 的 自由 能 


度 快 于 MM 的 1 


能 随时 间 的 变化 曲线 . 从 图 4 可 以 
在 体系 扩散 过 程 中 , 取向 差 角 较 大 的 相 界 与 
取向 差 角 较 小 的 相 界 扩散 体系 的 自由 能 均 呈 现 相 
同 的 下 降 变 化 趋势 , 这 符合 
自由 能 降低 的 要 求 . 但 是 对 比 图 
可 差 角 较 小 的 扩散 体系 
下 降 速 率 要 显著 快 于 MzMs 的 情形 , } 
能 相对 较 小 的 稳定 状态 . 对 于 取向 差 角 较 大 的 


里 闭 体系 互 扩散 过 程 中 
4a 和 b 可 以 看 出 , 取 
PF, Ms=Ms 情 况 下 , 自由 能 
日 首先 达到 


噬 速 


青 形 , 并 且 在 7=3x105At 左右 时 ， 


Ms=Ms 对 应 的 自由 能 达到 相对 稳定 状态 . MM 的 


自由 能 i 


已 


自由 能 的 下 降 速 率 随 Mwam 比 值 的 增 大 而 增 大 , 相 
对 于 取向 差 角 较 小 的 扩散 体系 ， 
定 状 态 时 的 时 间 缩 短 . 


到 相对 稳 


对 于 取向 差 角 较 小 的 扩散 体系 , 当 M=Ms 时 ， 
扩散 体系 不 存在 Kirkendall 效应, 不 存在 净 原 子 与 


禅 空位 及 由 


于 原子 空位 集聚 而 产生 的 


局 部 应 力 , 所 


金 属 


学 报 


以 在 封闭 的 扩散 体 


ChinaXiv 合 作 期 刊 


第 51 卷 


系 中 , 自由 能 由 相对 较 高 


值 逐 步 


降低 , 直至 达到 自由 能 相对 较 低 状态 . 而 当 MszMs 


时 , 由 于 Kirkendall 


效应 的 存在 , 扩散 体系 的 


局 部 区 


域 会 产生 一 定 的 六 


原子 与 净空 位 集聚 , 使 得 


域 出 现 应 力 , 抑制 


低 , 从 而 使 得 体系 的 自由 能 下 降 速率 要 小 于 
时 的 情形 . 此 外 , 由 于 扩散 体系 为 取向 差 角 


10000 


了 封闭 扩散 体系 的 自由 


能 的 降 


局 部 区 


Ma,=Ms 
相对 较 


ba 
9500 上 NA 
， 


9000 上 


Free energy 


8500 上 


8000 


2x105 3x105 4x105 5x105 6xl 
T 


Free energy 


| 
1x105 


图 4 不 同 原子 迁移 


率 下 , 取向 差 
的 扩散 过 程 体 系 的 自由 能 变化 


j 较 小 与 取向 差 


角 较 大 


Fig.4 Free energy evolution curves corresponding to differ- 


ent atomic mobilities for low misorientation (=2°) 


(a) and large misorientation (=12°) (b) (Ma and Ms 


are the atomic mobilities corresponding to A atom 


and B atom, respectively) 


10000 
N 
9500 | 


9000 上 上 


Free energy 


8500 上 


8000 上 


1x105 


图 5 M=100Ms 时 不 同 相 界 取向 差 角 扩 散 体系 中 自 


随时 间 的 变化 


2x105 3x105 4x105 5x10; 6x1l 
万 


05 


台 b 
El 


Fig.S Free energy evolution curves corresponding to dif- 


ferent misorientations for the same atomic mobility 


ratio (Ms=100Ms) 
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(a) (b) (0) 


fp 


图 6 模拟 结果 与 实验 结果 ”的 对 比 


EL 


Fig.6 Comparisons ofthe simulated results of Kirkendall void at the phase boundary with experiment observationp9 


(a) simulated result at low misorientation ((=2°) phase boundary 


(b) simulated result at large misorientation (=12°) phase boundary 


(c) experiment result™" 


小 的 扩散 体系 , 集聚 原子 与 空位 析出 与 潭 没 位 置 有 
限 , 所 以 随 着 A 和 B2 种 原子 迁移 率 差 值 的 增 大 , 体 
系 的 自由 能 下 降 速率 并 未 出 现 明显 的 差别 . 对 于 取 
向 差 角 较 大 的 扩散 体系 , 当 M=Ms 时 , 体系 的 自由 
能 下 降 原 因 类 似 于 取向 差 角 较 小 的 扩散 体系 的 情 
形 . 当 MMz#Ms 时 , 由 于 取向 差 角 较 大 的 相 界 处 空洞 
密度 增 大 , 集聚 的 原子 与 空位 可 以 在 更 多 位 置 析 出 
及 淹没, 从 而 随 着 A 和 B 2 种 原子 迁移 率 差 值 的 增 
大 , 由 Kirkendall 效 应 诱发 的 净空 位 与 净 原 子 扩散 
通 量 加 大 , 并 能 及 时 在 不 同 的 空洞 周围 发 生涯 没 及 
析出 , 从 而 使 得 自由 能 的 下 降 速 率 逐 渐 增 大 . 

图 5 反映 的 是 M=100Ms 时 ,不同 取 向 差 角 扩散 
体系 中 自由 能 随时 间 的 变化 曲线 . 可 以 看 出 , 不 同 
取向 差 角 的 扩散 体系 中 , 自由 能 均 由 相对 较 大 值 不 
断 降 低 . 随 着 相 界 取向 夹 角 的 增 大 , 体系 的 自由 能 
下 降 速 率 逐 渐 增 大 . 这 是 因为 相 界 取向 差 角 增 大 ， 
相 界 处 位 错 核 (在 本 研究 中 可 以 作为 空位 的 淹没 源 ) 
的 密度 增 大 , 在 随后 的 扩散 过 程 中 形成 的 空洞 的 密 
度 也 逐渐 增 大 , 集聚 的 原子 与 空位 可 以 在 更 多 的 空 
洞 周 围 析出 及 淹没 , 集聚 的 应 力 迅 速 得 到 释放 , 从 
而 使 得 体系 的 自由 能 下 降 速 率 加 快 . 在 自由 能 演化 
的 后 期 , 随 着 空洞 的 纵向 移动 及 横向 连通 , 体系 的 
自由 能 进一步 降低 , 直至 达到 自由 能 相对 较 低 的 最 
终 稳 定 状态 . 

3 与 实验 的 对 比 

图 6a 和 b 分 别 是 取向 差 角 较 小 及 取向 差 角 较 大 

的 相 界 空洞 的 模拟 结果 , 图 6c 是 Cu/NiAl 扩 散 侦 在 


850 CC 退火 保温 54h 后 实际 观察 到 的 相 界 空洞 中 
对 比 图 6a 与 c 中 数字 1 所 示 区 域 , 可 以 验证 模拟 得 
到 的 单个 孤立 空洞 的 准确 性 , 对 比 图 6b 与 c 中 数字 
2 所 示 区 域 , 可 以 验证 模拟 空洞 连通 的 准确 性 . 上 述 
对 比 证 明了 模拟 结果 真实 可 靠 . 

4 结论 

(1) 对 于 取向 差 角 较 小 的 相 界 (9=2s), 相 界 处 
形成 的 空洞 密度 较 低 , 各 空洞 之 间 影 响 较 小 ; 对 于 
取向 差 角 较 大 的 相 界 (9=12°), 相 界 处 空洞 密度 较 
大 , 各 空洞 之 间 影 响 较 大 . 单个 微 空洞 在 演化 过 程 
中 向 高 原子 迁移 率 一 侧 (a 相 ) 移动 , 空洞 的 形状 由 
最 初 的 平行 四 边 形 向 六 角形 演化 . 

(2) 取向 差 角 较 小 的 相 界 空洞 扩展 过 程 伴 有 相 
界 的 移动 (向 高 原子 迁移 率 一 侧 (w 相 ) 移动 ), 使 B 相 
逐渐 长 大 , w 相 逐渐 缩小 . 取向 差 角 较 大 的 相 界 空洞 
扩展 过 程 中 还 会 沿 相 界 方向 扩展 , 使 空洞 连通 , 将 
相 界 割裂 开 , 割裂 的 两 侧 相 界 呈现 锯齿 状 . 

(3) 在 封闭 的 互 扩散 体系 中 , 体系 的 自由 能 由 相 
对 较 大 值 逐 渐 降 低 . 对 于 取向 差 角 较 小 的 相 界 , 随 
着 MMs 比值 的 增 大 , 体系 的 自由 能 下 降 速 率 没有 
明显 的 差异 . 对 于 取向 差 角 较 大 的 相 界 , 随 着 MMs 
比值 的 增 大 , 互 扩散 体系 的 自由 能 下 降 速率 逐渐 增 
大 . 随 着 取向 夹 角 的 增 大 , 演化 过 程 中 空洞 的 密度 
增 大 , 集聚 的 空位 与 原子 在 更 多 位 置 发 生 淹没 , 自 
由 能 的 下 降 速率 逐渐 增 大 . 


Ll 
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